Entanglement entropy in all dimensions 
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It has long been conjectured that the entropy of quantum fields across boundaries scales as the 
boundary area. This conjecture has not been easy to test in spacetime dimensions greater than four 
because of divergences in the von Neumann entropy. Here we show that the Renyi entropy provides 
a convergent alternative, yielding a quantitative measure of entanglement between quantum field 
theoretic degrees of freedom inside and outside hypersurfaces. For the first time, we show that the 
entanglement entropy in higher dimensions is proportional to the higher dimensional area. We also 
show that the Renyi entropy diverges at specific values of the Renyi parameter q in each dimension, 
but this divergence can be tamed by introducing a mass of the quantum field. 
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Entanglement, a term first coined by Schrddinger, is 
an intriguing and quintessentially quantum mechanical 
property, which correlates microscopic systems in a pre- 
cise way, even if they are separated by large distances. On 
the one hand it gives rise to apparent contradictions (such 
as the EPR paradox) and on the other, hides enormous 
untapped resources for computation and communication 
(e.g., teleportation). A mathematically precise way of 
measuring entanglement has remained elusive however, 
except in the simplest cases where the combined system 
is in a pure state, i.e., for which all quantum numbers 
are known. Usually, the entanglement entropy is com- 
puted as the von Neumann entropy associated with p: 
5 vN = — tr(plnp). For recent reviews, see [1, 2]. 

Entanglement as a source of entropy has also been ex- 
amined in an entirely different context, that of the mi- 
croscopic origin of black hole entropy. The basic idea is 
that the long-range entanglement of the quantum fields 
across a black-hole horizon can leave its mark on the re- 
duced density matrix of the external degrees of freedom 
(DOF) which, in turn, accounts for the black-hole en- 
tropy. Starting with [3, 4], it was demonstrated that to 
leading order, and for a scalar field in its global ground 
state, the entanglement entropy between its DOF inside 
and outside a black hole horizon is proportional to the 
area of the horizon. Recently, it was demonstrated that 
the DOF at or near the horizon contribute to most of 
this entropy [5] and that excited states of the field lead 
to sub- leading order power law corrections [6, 7]. 

However, entanglement as a source of black-hole en- 
tropy has a couple of drawbacks: (i) The proportionality 
constant depends on the ultra-violet cut-off and the num- 
ber of fields present (also these are in general independent 
of each other, although it was recently heuristically ar- 
gued that the requirement of the stability of the cosmos 
related the two and naturally related the cut-off to the 
Planck mass [8]), and (ii) it is not evident that the area 



proportionality holds for higher spacetime dimensions, 
say D + 2 > 4 (throughout D + 1 denotes the number of 
spatial dimensions)? 

Although there have been attempts to obtain uni- 
versal expressions for higher dimensions from the two- 
dimensional entropy c-functions [1], such is yet to be 
shown from first principles. For higher dimensional 
black-holes, the brick-wall entropy contains extra diver- 
gent terms other than from the ultra-violet modes [9]. 
Following Srcdnicki [4], if one regularizes the entropy 
function by introducing a radial lattice, the sum of par- 
tial waves does not converge and the entropy turns out to 
be infinite in higher dimensions. Here, we show that the 
entropy-area relation can be obtained for all dimensions 
by using a different measure of entropy, and furthermore, 
that the divergences are similar in nature to the infrared 
divergences in QED, which can be tamed by introducing 
a mass to the field. 

While the von Neumann entropy is the most common 
measure of entanglement, it is neither the most general, 
nor unique. There are other measures, such as the Renyi 
and Tsallis entropies, which under certain limits reduce 
to the von Neumann entropy. In this work we study 
entanglement via the Renyi entropy defined as 



(i) 



In the limit that the parameter q — > 1, the Renyi en- 
tropy reduces to its von Neumann counterpart. Also, 
like von Neumann entropy (and unlike Tsallis entropy) 
Renyi entropy is additive and has maximum value ln(n) 
for pi = 1/n. 

Consider a free, massless real scalar field propagating 
in a (D + 2)-dimensional flat spacetime with action 

~ J dtdr r D dn D (r] ab d a <f> d b <P + g 6 " 6 " d e J> d e J>), (2) 
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where the Latin indices (a and b) take the values t and r, 
flu is the D-dimcnsional solid angle, rf> is the Minkowski 
metric, and are the metric coefficients for the an- 

gular coordinates. Decomposing the scalar field in terms 
of real hyper-spherical harmonics one has 
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where i G {1, . . . , D — 1}. The canonical momentum 
of the field is then ire mi — d t (pe mi and the Hamiltonian 

H = Y,i, mi H t,mi where 
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Note that the Hamiltonian of a massless scalar field prop- 
agating in a general (i.e., not necessarily flat) spherically 
symmetric spacetime at fixed Lemaitre time coordinate 
reduces to the above Hamiltonian [7] . 

For simplicity, we assume that the scalar field is in the 
ground (vacuum) state. As mentioned earlier, our inter- 
est is in determining the quantity of information shared 
by modes across a "horizon." This may be found by elim- 
inating (tracing out) the quantum degrees of freedom as- 
sociated with the scalar field inside a spherical region of 
radius 1Z (the location of our "horizon"). The resulting 
reduced density matrix can then be used to determine 
the strength of quantum correlations across the horizon. 

It is not possible to obtain a closed form analytic ex- 
pression for the density matrix and hence, we need to 
resort to numerical methods. In order to do that we 
take a spatially uniform radial grid, {rj}, with a = 
rj+i — rj. To achieve better precision, we adopt the 
middle-prescription in discretizing the terms containing 
the derivatives, so 

A(r)d r [G{r)] A(j + \) G{j ± ^ ° {]) , (5) 

Discretizing the Hamiltonian (4), and suppressing the 
subscripts l,m,i, leads to 

1 N 

H = ^H ] = — 5^(5y7r? + if j Kij (p^, (6) 

3 i,3 

where the interaction matrix, , is given by 
[{] + \)° +f- 2 m + D-l)]5)8] 
+ jo [U + l f+ U - \) D + 3 D - 2 l(l + D - 1)] Jjtf 

+ I) fij+l (' + I) xi+1 
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and 2 < k < (N — 1). The procedure to obtain the 
entanglement entropy in higher dimensions is similar to 
the one discussed in Refs. [4, 10]. In this work we assume 
that the quantum state corresponding to the Hamiltonian 
of the iV-Harmonic oscillator system (6) is the ground 
state with wave-function &gs(xi, ■ ■ ■ ,x n ;t\, . . . ,ijv_„). 
The reduced density matrix p(t, t ) is obtained by tracing 
over the first n of the N oscillators 
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(ng=i da*) *gs(^i » • • • » x n ;t)^Gs{xi,. ..,x n ; t ). (8) 



The Renyi entanglement entropy in (D + 2)-dimensional 
spacetime is then given by 

^(n.JVj^fSf + fl-ljWW^^JV), (9) 



where (n, N) is the Renyi entropy for partial waves 
with total angular momentum I and the angular degen- 
eracy factor W(£) is given by [9, 11] 



(£ + D-2)l 
{D-l)\£\ 
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To highlight the advantage of using the Renyi entropy as 
the measure of the entanglement entropy as compared to 
von Neumann entropy, we determine the Renyi entropy 
in the large £ 3> N limit. In this limit, Eq. (9) becomes 



S(«)(n,JV)~2(2^ + D-l)^e«, (11) 



where the asymptotic eigenvectors and angular degener- 
acy factors, respectively, are 
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Substituting these expressions into Eq. (11) yields 

S iq \ ni N) ~^£ D ~ 1 - 4q . (13) 

e 

In the limit Eq. (13) is identical to the asymptotic 

limit of the von Neumann entropy, and diverges for all 
D > 2 (also as noted in [4]). By contrast, if that limit is 
not taken, then Eq. (13) converges for all q > D/4. For 
instance, in 4-dimensional spacetime this implies that for 
all values of q > \ entanglement entropy converges, while 
in 6-dimensions, this entropy converges only for q > 1. 
Herein lies the main advantage of using Renyi entropy. 
The above asymptotic expression also provides an under- 
standing as to why von Neumann entropy converges for 
4-dimcnsions and not for 6-dimensions. 

We will now show that the Renyi entropy indeed pro- 
vides a good measure for the entanglement entropy in all 
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dimensions for q > D/4, and also reproduces the area 
proportionality in each case. We compute the Renyi en- 
tropy numerically for the discretizcd Hamiltonian (6). 
The computations are done using Matlab for the lattice 
size N — 300, 100 < n < 200 and the relative error in the 
computation of Renyi entropy is 10~ 6 . The computations 
were done for 4-, 5-, 6- and 10-dimensional spacetimes. 
However, in this manuscript we only present the results 
for 5 and 10 dimensions. 
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FIG. 2: Log-log plot of Renyi entropy, S^ q \ versus the scaled 
radius of the sphere, s = TZ/a = n + =, for 10-dimensional 
spacetime. The lattice size is N = 300 and 100 < n < 200. 
The fits are: S (2 ' 5) = 6.08 x 10 9 s 8 ' 00 , S w = 5.71 x 10 9 s 8 ' 00 , 
S (5) = 5.35 x 10 9 s 8 - 00 , S (6) = 5.14 x 10 9 s 8 - 00 , S (6 ' 5) = 5.14 x 
10 9 s 8.o 0j 5 (7) =5 .oixl0 9 s 8 00 ,S (8) = 4.89xl0 9 s 8 - 00 , S (9) = 
4.81 x 10 9 s 8 00 , 5 (9 5) = 4.60 x 10 9 s 8 - 00 , and S (10) = 4.76 x 
10 9 s 8 ' 00 . 
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FIG. 1: Log-log plot of Renyi entropy, S^ q \ versus the 
scaled radius of the sphere, s = TZ/a — n+i, for 5- 
dimensional spacetime. The lattice size is N — 300 and 
100 < n < 200. The black dots represent the numerical 
output and straight lines denote lines of best fit. These fits 
are: 5 (0 9) = 0.69s 2 ", = 0.26s 2 ", S (2) = 0.068s 2 ", 



S< 3 > = 0.056 s 2 ' 98 , SW = 
S (10) = 0.047s 2 - 98 , 5 (20) 
and S (100) = 0.045 s 2 - 97 . 
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In Figs. 1 and 2 we have plotted lnS^ versus \n(R/a) 
for 5- and 10-dimensional spacetimes, respectively. From 
the best fit curves in Fig. 1 we see that the Renyi entropy 
for different values of q generically scales with an approx- 
imate power law ~ (TZ/a) 3 . Although as we increase 
<7, the prefactor decreases, nonetheless the power remains 
close to 3. Similarly, in the case of 10-dimensional space- 
time (Fig. 2), for all q, the power always remains close 
to 8. It is important to note that, although it has long 
been conjectured in the literature that the entanglement 
entropy is proportional to the area of the D-dimcnsional 
surface (for D > 3), to our knowledge this is the first 
time this has been explicitly shown in higher spatial di- 
mensions. 

In Fig. 3 the prefactor S'profactor, (the ratio of the 
Renyi entropy to 5-dimensional area) is plotted versus 
q, from which we infer the following: (i) Renyi entropy 
saturates to a constant value for large q. (ii) Consistent 
with the asymptotic analysis of Eq. (13), the Renyi en- 
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FIG. 3: Plot of the prefactor, Sprcfactor, of the Renyi entropy 
in 5-dimensional spacetime versus \n(q) . The small circles rep- 
resent the numerical output for each run; the red line shows 
a spline interpolation. The divergence is apparent as q — > |. 



tropy diverges logarithmically as q — > |. Similarly, for 
(D + 2)-dimensional spacetime, it diverges as q — > D/4. 

To understand this divergence further, we now consider 
a scalar field of mass m. The new interaction matrix is 
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where Kij is given by Eq. (7). In Fig. 4 we plot In S^ q ' vs 
ln(i?/et) again for 5-dimensional spacetime but for q = |, 
and for different values of to, and show that the diver- 
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FIG. 4: Log-log plot of Renyi entropy, 5 , ' 9 \ versus the scaled 
radius of the sphere, s = IZ/a = n + i, in 5-dimensional 
spacetime for q = | (which is divergent in the massless case) 
and different non-zero masses (er = ma. The lattice size is 
AT = 300 and 100 < n < 200. The black dots represent the 
numerical output and straight lines are the lines of best fit. 
These fits are: S (3/4) (m = 1) = 4.25 s 2 ' 97 , S (3/4) (m = 1.5) = 
4.16 s 2 ", SW A \m = 10) = 3.20s 2 - 99 and S< 3 / 4 >(m = 50) = 
2.07 s 2 ' 99 . 



ularly disturbing since it is a free field theory without 
interactions, for which one normally understands how to 
absorb those divergences. Here on the other hand, the 
ultra-violet cut-off seems real, without however any clear 
physical meaning or knowledge about how it is imple- 
mented in the real world. 

3. Over the last decade it has been shown that higher- 
dimensional black-holes can have a much richer topolog- 
ical structure than 4-dimensional black-holes. For in- 
stance, it was shown by Emparan and Reall [12] that 
horizons with a topology of S 2 x S 1 can exist for some 
asymptotically flat spacetimes in 5-dimensions. It would 
be interesting to investigate whether the entanglement 
entropy-area relation holds for such topologies. 
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gences disappears! We also note that as the mass in- 
creases, the prefactor S'pretactor of the Renyi entropy mod- 
estly decreases. 

To summarize, in this work we have computed 
Renyi entropy after tracing out the degrees of freedom 
inside a D-dimensional sphere of radius 1Z, and show that 
the entanglement entropy is proportional to the area of 
the D-dimensional surface (for D > 3). Although conjec- 
tured before, this is the first time to our knowledge, that 
such area proportionality has been shown explicitly for 
higher spatial dimensions. We have also shown that the 
logarithmic divergence of Renyi entropy can be tamed by 
considering massive fields. 

We conclude with the following points: 1. Eq. (9) 
involves summation over partial waves £ that extend to 
infinity. Since the massless Klein-Gordon field does not 
have an intrinsic length scale, all angular momentum 
modes may contribute. However, for large enough values 
of £, the energy corresponding to the angular momentum 
modes will exceed the Planck energy. One way to address 
this difficulty may be to introduce higher-order deriva- 
tives to the scalar field action. However, it is not clear 
whether the entanglement entropy-area relation would 
emerge in this case. 

2. The entanglement entropy of a massless scalar- field 
is plagued by ultra-violet divergences. This is partic- 
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